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We propose a comprehensive theory of dark matter that explains the recent proliferation of un- 
expected observations in high-energy astrophysics. Cosmic ray spectra from ATIC and PAMELA 
require a WIMP with mass M x ~ 500 — 800 GeV that annihilates into leptons at a level well above 
that expected from a thermal relic. Signals from WMAP and EGRET reinforce this interpretation. 
Limits on p and 7r°-7's constrain the hadronic channels allowed for dark matter. Taken together, we 
argue these facts imply the presence of a new force in the dark sector, with a Compton wavelength 
m^ 1 ^ 1 GeV -1 . The long range allows a Sommerfeld enhancement to boost the annihilation cross 
section as required, without altering the weak-scale annihilation cross section during dark matter 
freeze-out in the early universe. If the dark matter annihilates into the new force carrier cj>, its low 
mass can make hadronic modes kinematically inaccessible, forcing decays dominantly into leptons. 
If the force carrier is a non-Abelian gauge boson, the dark matter is part of a multiplet of states, and 
splittings between these states are naturally generated with size am^ ~ MeV, leading to the eX- 
citing dark matter (XDM) scenario previously proposed to explain the positron annihilation in the 
galactic center observed by the INTEGRAL satellite; the light boson invoked by XDM to mediate 
a large inelastic scattering cross section is identified with the <f> here. Somewhat smaller splittings 
would also be expected, providing a natural source for the parameters of the inelastic dark matter 
(iDM) explanation for the DAMA annual modulation signal. Since the Sommerfeld enhancement 
is most significant at low velocities, early dark matter halos at redshift ~ 10 potentially produce 
observable effects on the ionization history of the universe. Because of the enhanced cross section, 
detection of substructure is more probable than with a conventional WIMP. Moreover, the low veloc- 
ity dispersion of dwarf galaxies and Milky Way subhalos can increase the substructure annihilation 
signal by an additional order of magnitude or more. 

PACS numbers: 95.35.+d 



Thermal WIMPs remain one of the most attractive candidates for dark matter. In addition to appearing generically 
in theories of weak-scale physics beyond the standard model, they naturally give the appropriate relic abundance. Such 
particles also are very promising in terms of direct and indirect detection, because they must have some connection 
to standard model particles. 

Indirect detection is particularly attractive in this respect. If dark matter annihilates to some set of standard model 
states, cosmic ray detectors such as PAMELA, ATIC and Fermi/GLAST have the prospect of detecting it. This is 
appealing, because it directly ties the observable to the processes that determine the relic abundance. 

For a weak-scale thermal particle, the relic abundance in the case of s-wave annihilation is approximately set by 



For perturbative annihilations, s-wave dominates in the late universe, so this provides an approximate upper limit on 
the signal that can be observed in the present day. Such a low cross section makes indirect detection, whereby the 
annihilation products of dark matter are detected in cosmic ray detectors, a daunting task. 



I. PAMELA/ ATIC AND NEW DARK FORCES 
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However, recent experiments have confirmed the longstanding suspicion that there are more positrons and electrons 
at lOs-lOOs of GeV than can be explained by supernova shocks and interactions of cosmic ray protons with the ISM. 
The experiments are 

• PAMELA 

The Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics has reported results [l| indicating 
a sharp upturn in the positron fraction (e + /(e + + e - )) from 10 — 100 GeV, counter to what is expected from 
high-energy cosmic rays interacting with the interstellar medium (ISM). This result confirms excesses seen in 
previous experiments, such as HEAT 0, 0| and AMS-01 Q]. One possible explanation for this is dark matter 
annihilation into e + e , but this requires a large cross section [8| . 

• ATIC 

The Advanced Thin Ionization Calorimeter is a balloon-borne cosmic ray detector which studies electrons 
and positrons (as well as other cosmic rays) up to ~ TeV energies, but cannot distinguish positrons and 
electrons. The primary astrophysical sources of high-energy electrons are expected to be supernovae: electrons 
are accelerated to relativistic speeds in supernova remnants and then diffuse outward. The ATIC-2 experiment 
reported a 4 — 6a excess (over a simple power law) in its e + + e~ data Q at energies of ~ 300 — 800 GeV, with 
a sharp cutoff in the 600 — 800 GeV range. Dark matter would seem a natural candidate for this as well, with 
its mass scale determining the cutoff. 

• WMAP 

Studies of the WMAP microwave emission from the galactic center show a hard component not spatially corre- 
lated with any known galactic emission mechanism. This "WMAP haze" 0,0] can be explained as synchrotron 



radiation from electrons and positrons produced from dark matter annihilation in the galactic center 12 1 



• EGRET 



Gamma-ray measurements in the galactic center (inner 5°) provide hints of an excess at 10 — 50 GeV 1_3| . Strong 



et al. reanalyzed EGRET data and found a harder spectrum at these energies than previously derived, using 



the improved EGRET sensitivity estimates of 14[. Despite poor spatial resolution, Strong ct al. found an 



excess in this energy range above the expected ir° gamma-ray emission from cosmic ray protons interacting with 



the interstellar medium (see Fig. 8 of 13(). Such 7's could naturally arise from inverse-Compton scattering of 



high-energy electrons and positrons off of starlight and the cosmic microwave background (CMB). 

Taken together, these make a compelling case for excessive electronic production in the galaxy. While individual 
astrophysical explanations may exist for each signal (pulsar wind nebulae for PAMELA and ATIC 15, 1^, for 
instance, supernovae for the WMAP haze |63|), the data cry out for a unified explanation. Dark matter annihilations 
provide an appealing candidate. 

In addition to the above, there are two other anomalies that are worth mentioning. The INTEGRAL 511 keV 
signal indicates ~ 3 x 10 42 e + /s annihilating in the galactic center, far more than expected from supernovae. The 



spectrum suggests that these positrons are injected with relatively low energies (E ^ few MeV) , and so form a distinct 



population from those above. eXciting dark matter (XDM) [l8j can naturally explain such low energy positrons with 
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~ 1 MeV excited states of the dark matter, while still producing the high-energy positrons via annihilation 
Lastly, there is the DAM A/LIBRA indication of an annual modulation consistent with that expected from dark matter 
induced nuclear scattering [2l|. Such a signal is difficult to reconcile with null results of other exp eriments, but can 
be reconciled with ~ 100 keV excited states in the "inelastic dark matter" (iDM) scenario 

HOB. Although we 



are motivated by the specific signals above (PAMELA/ ATIC as well as the haze and EGRET), the picture we are 
led to for explaining them naturally incorporates the necessary ingredients to explain the INTEGRAL and DAMA 
signals as well. 

Focusing on only the high-energy positrons and electrons, there are a number of challenges to any model of dark 
matter. These are: 



• A large cross section 

Studies of PAMELA and ATIC signals seem to require a cross section much larger than what is allowed by 
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thermal relic abundance. Boost factors of 0(100) or more above what would be expected for a thermal WIMP 
are required to explain these excesses 
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• A large cross section into leptons 

Typical annihilations via Z bosons produce very few hard leptons. Annihilations into W bosons produce hard 
leptons, but many more soft leptons through the hadronic shower. Higgs bosons and heavy quarks produce even 
softer spectra of leptons, all of which seem to give poor fits to the data. At the same time, absent a leptophilic 
gauge boson, it is a challenge to construct means by which dark matter would annihilate directly to leptons. 

• A low cross section into hadrons 

Even if a suitably high annihilation rate into leptons can be achieved, the annihilation rate into hadronic modes 
must be low. Limits from diffuse galactic gamma rays [27| . as well as gamma rays from the center of the 
galaxy constrain the production of 7r°'s arising from the dark matter annihilation. PAMELA measurements of 
antiprotons tightly constrain hadronic annihilations as well . Consequently, although quark and gauge boson 
annihilation channels may occur at some level, the dominant source of leptons must arise through some other 
channel. 

The combination of these issues makes the observed high-energy anomalies - especially ATIC and PAMELA 
difficult to explain with thermal dark matter annihilation. However, we shall see that the inclusion of a new force in 
the dark sector simultaneously addresses all of these concerns. 

New Forces in the Dark Sector 

A new interaction for the dark sector can arise naturally in a variety of theories of physics beyond the standard 
model, and is thus well motivated from a theoretical point of view. Although there are strong limits on the self- 



interaction scattering cross section from structure formation 28|, |29( , the presence of some new force-carrying boson 



should be expected, with only the mass scale in question. A light boson could arise naturally if its mass scale is 



generated radiatively [30| , or if it were a pseudo-goldstone boson. 

One of the important modifications that can arise with a new light boson is an enhancement of the annihilation 
cross section via a mechanism first described by Sommerfeld [3~J ] . The presence of a new force carrier can distort the 
wave function of the incoming particles away from the plane-wave approximation, yielding significant enhancements 



(or suppressions) to annihilation cross sections 6J|. Equivalently, ladder diagrams involving multiple exchanges of 
the force carrier must be resummed (Fig. 1). As we shall describe, the Sommerfeld enhancement can only arise if the 
gauge boson has a mass m^^S aMjjM ~ few GeV. Thus, with the mass scale of ~ 800 GeV selected by ATIC, and 
the large cross sections needed by both ATIC and PAMELA, the mass scale for a new force carrier is automatically 
selected. Interactions involving W and Z bosons are insufficient at this mass scale. 

Once this new force carrier <fi is included, the possibility of a new annihilation channel XX ~~ > 00 opens up, which can 
easily be the dominant annihilation channel. Absent couplings to the standard model, some of these particles could 
naturally be stable for kincmatical reasons; even small interactions with the standard model can then lead them to 
decay only into standard model states. If they decay dominantly into leptons, then a hard spectrum of positrons arises 
very naturally. Motivated by the setup of XDM, Cholis, Goodenough and Weiner [3| first invoked this mechanism of 
annihilations into light bosons to provide a simple explanation for the excesses of cosmic ray positrons seen by HEAT 
without excessive antiprotons or photons. Simple kinematics can forbid a decay into heavier hadronic states, and as 
we shall see, scalars lighter than ~ 250 MeV and vectors lighter than ~ GeV both provide a mode by which the dark 
matter can dominantly annihilate into very hard leptons, with few or no 7r°'s or antiprotons. 

In the following sections, we shall make these points more concretely and delineate which ranges of parameters most 
easily explain the data; but the essential point is very simple: if the dark matter is O(800 GeV) and interacts with 
itself via a force carrier with mass ~ GeV, annihilation cross sections can be considerably enhanced at present 
times via a Sommerfeld enhancement, far exceeding the thermal freeze-out cross section. If that boson has a small 
mixing with the standard model, its mass scale can make it kinematically incapable of decaying via a hadronic shower, 
preferring muons, electrons and, in some cases, charged pions, and avoiding constraints from 7r°'s and antiprotons. 
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If the force-carriers are non-Abelian gauge bosons, we shall see that other anomalies may be incorporated naturally 
in this framework, explaining the INTEGRAL 511 keV line via the mechanism of "eXciting dark matter" (XDM), 
and the DAMA annual modulation signal via the mechanism of "inelastic dark matter" (iDM) . We shall see that the 
excited states needed for both of these mechanisms [0(1 MeV) for XDM and 0(100 keV) for iDM] arise naturally with 
the relevant mass splitting generated radiatively at the correct scale. 



II. SOMMERFELD ENHANCEMENTS FROM NEW FORCES 

A new force in the dark sector can give rise to the large annihilation cross sections required to explain recent data, 



through the "Sommcrfcld enhancement" that increases the cross section at low velocities [65j. A simple classical 
analogy can be used to illustrate the effect. Consider a point particle impinging on a star of radius R. Neglecting 
gravity, the cross section for the particle to hit and be absorbed by the star is <7o = irR 2 . However, because of gravity, 
a point coming from a larger impact parameter than R will be sucked into the star. The cross section is actually 
a = Tr^as' where b max is the largest the impact parameter can be so that the distance of closest approach of the 
orbit is R. If the velocity of the particle at infinity is v, we can determine b max trivially using conservation of energy 
and angular momentum, and we find that 



where v 2 sc = 2G^M/R is the escape velocity from the surface of the star. For v <C v esc , there is a large enhancement 
of the cross section due to gravity; even though the correction vanishes as gravity shuts off (Gn — > 0), the expansion 
parameter is 2GnM/(Rv 2 ) which can become large at small velocity. 

The Sommerfeld enhancement is a quantum counterpart to this classical phenomenon. It arises whenever a particle 
has an attractive force carrier with a Compton wavelength longer than {olMdm) , i-e. dark matter bound states 
are present in the spectrum of the theory. (We generically refer to a ~ coupling 2 /47r, assuming such couplings are 
comparable to those in the standard model, with 10~ 3 ^5 a'S 10 _1 .) 

Let us study this enhancement more quantitatively. We begin with the simplest case of interest, namely, a particle 
interacting via a Yukawa potential. We assume a dark matter particle x coupling to a mediator <f> with coupling 
strength A. For s-wave annihilation in the nonrelativistic limit, the reduced two-body wavefunction obeys the radial 
Schrodinger equation, 

—ip"(r) - V(rWr) = -m y v 2 rp(r), (3) 

where the s-wave wavefunction ^(r) is related to ■0( r ) a s ^(r) = ip{r)/r, v is the velocity of each particle in the 
center-of-mass frame (here we use units where h = c = 1), and for scalar 4> the potential takes the usual Yukawa form, 

V(r) = --^- e - m * r . (4) 

The interaction in the absence of the potential is pointlike. As reviewed in the appendix, the Sommerfeld enhancement 
in the scattering cross section due to the potential is given by 



S = 



(5) 



where we solve the Schrodinger equation with boundary conditions ip(0) = 0, ip{r) — > sin(fcr + S) as r — » oo. In the 
recent dark matter literature, a different but completely equivalent expression is used, with 

SH<Moo)M0)| 2 (6) 

where we solve the Schrodinger equation with the outgoing boundary condition ^'(oo) = im x vip(oo) (3^ |. 
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Defining the dimensionless parameters 




(7) 



(8) 



In the limit where the <fi mass goes to zero (e^ — > 0), the effective potential is just the Coulomb potential and Eq. 
[8] can be solved analytically, yielding an enhancement factor of, 



For nonzero m$ and hence nonzero e^, there are two important qualitative differences. The first is that the Sommerfeld 
enhancement saturates at low velocity-the attractive force has a finite range, and this limits how big the enhancement 
can get. At low velocities, once the deBroglie wavelength of the particle (Mv)' 1 gets larger than the range of the 
interaction ml , or equivalently once e v drops beneath e^, the Sommerfeld enhancement saturates at S ~ The 
second effect is that for specific values of e^, the Yukawa potential develops threshold bound states, and these give rise 
to resonant enhancements of the Sommerfeld enhancement. We describe some of the parametrics for these effects in 
the appendix, but for reliable numbers Eq. [8] must be solved numerically, and plots for the enhancement as a function 
of £0 and e v are given there. As we will see in the following, we will be interested in a range of ~ 100 MeV - GeV; 
with reasonable values of a, this corresponds to in the range ~ 10~ 2 — 10 _1 , yielding Sommerfeld enhancements 
ranging up to a factor ~ 10 3 — 10 4 . At low velocities, the finite range of the Yukawa interaction causes the Sommerfeld 
enhancement to saturate, so the enhancement factor cannot greatly exceed this value even at arbitrarily low velocities. 
The nonzero mass of the (j) thus prevents the catastrophic overproduction of gammas in the early universe pointed 



Having obtained the enhancement S as a function of e v and e^, we must integrate over the velocity distribution 
of the dark matter in Earth's neighborhood, to obtain the total enhancement to the annihilation cross section for 
a particular choice of <f> mass and coupling A. We assume a Maxwell-Boltzmann distribution for the one-particle 
velocity, truncated at the escape velocity: 



The truncation does not significantly affect the results, as the enhancement factor drops rapidly with increasing 
velocity. The one-particle rms velocity is taken to be 150 km/s in the baseline case, following simulations by Governato 
et al. [34j | . Fig. [T] shows the total enhancement as a function of m l p/m x and the coupling A for this case. 

One can inquire as to whether dark matter annihilations in the early universe experience the same Sommerfeld 
enhancement as dark matter annihilations in the galactic halo at the present time. This is important, because we are 
relying on this effect precisely to provide us with an annihilation cross section in the present day much larger than 
that in the early universe. However, it turns out that particles leave thermal equilibrium long before the Sommerfeld 
enhancement turns on. This is because the Sommerfeld enhancement occurs when the expansion parameter a/v = 1 / e v 
is large. In the early universe, the dark matter typically decouples at Tqmb ~ TO x/^0, or u ~ 0.3c. Since we are 
taking a^S 0.1 generally, the Sommerfeld enhancement has not turned on yet. More precisely, in the early-universe 
regime e v ^> e^, so we can use Eq. [H] for the massless limit. Where e v 3> 1, Eq. [5] yields S« 1 + 7r/2e„: thus the 
enhancement should be small and independent of m^. Fig. [2] shows this explicitly. We are left with the perturbative 
annihilation cross section a ~ a 2 /m 2 which gives us the usual successful thermal relic abundance. 

At some later time, as the dark matter velocities redshift to lower values, the Sommerfeld enhancement turns on 
and the annihilations begin to scale as a -5 / 2 (before kinetic decoupling) or a~ 2 (after decoupling). From decoupling 
until matter-radiation equality, where the Hubble scale begins to evolve differently, or until the Sommerfeld effect is 
saturated, dark matter annihilation will produce a uniform amount of energy per comoving volume per Hubble time. 



S= |^(oo)/^(0)| 2 



ir/e v 



(9) 
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out by [33 1 . 




(10) 
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FIG. 1: Contours for the Sommerfeld enhancement factor S as a function of the mass ratio m,j > /m x and the coupling constant 
\, a — 150 km/s. 



This uniform spread of energy injected could have potentially interesting signals for observations of the early universe. 
An obvious example would be a possible effect on the polarization of the CMB, as described in 35|, |3y, l37f. Because 
at the time of matter-radiation equality, the dark matter may have slowed to velocities of v ~ 10 _6 c or slower, the 
large cross sections could yield a promising signal for upcoming CMB polarization observations, including Planck. 
However, we emphasize that saturation of the cross section at low v avoids the runaway annihilations discussed by 

El. 



III. MODELS OF THE SOMMERFELD FORCE AND NEW ANNIHILATION CHANNELS 

What sorts of forces could give rise to a large Sommerfeld enhancement of the dark matter annihilation? As we 
have already discussed, we must have a light force carrier. On the other hand, a massless particle is disfavored by the 
agreement between big bang nucleosynthesis calculations and primordial light element measurements (38| . as well as 
constraints from WMAP on new relativistic degrees of freedom j^. Thus, we must have massive degrees of freedom, 
which are naturally light, while still coupling significantly to the dark matter. There are three basic candidates. 

• The simplest possibility is coupling to a light scalar field, which does give rise to an attractive interaction. 
However, given that we need an 0(1) coupling to the DM fields, this will typically make it unnatural for the 
scalar to stay as light as is needed to maximize the Sommerfeld enhancement, unless the dark matter sector 
is very supersymmetric. This can be a challenge given that we are expecting the dark matter to have a mass 
O(TeV). Consequently, the natural scale for a scalar which couples to it would also be O(TeV), although this 
conclusion can be evaded with some simple model-building. 

• The scalar could be naturally light if it is a pseudoscalar ir with a goldstone-like derivative coupling to matter 
1 / FJ^d^n. This does lead to a long-range spin-dependent potential of the form V[r) = i {Si • §2 — 3 Si • f§2 ■ r) , 
but the numerator vanishes when averaged over angles, so there is no long-range interaction in the s-wave and 
hence no Sommerfeld enhancement. 
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FIG. 2: Contours for the Sommerfeld enhancement factor S as a function of the mass ratio m,j > /m x and the coupling constant 
A, at a temperature Tcmb = m x /20 (the enhancement is integrated over a Boltzmann distribution with a = 0.3c). 



• Finally, we can have a coupling to spin-1 gauge fields arising from some dark gauge symmetry Gdark- Since 
the gauge fields must have a mass O(GeV) or less, one might worry that this simply begs the question, as the 
usual explanation of such a light gauge boson requires the existence of a scalar with a mass of O(GeV) or less. 
However, because that scalar need not couple directly to the dark matter, it is sufficiently sequestered that its 
small mass is technically natural. Indeed, the most straightforward embedding of this scenario within SUSY 
[jot naturally predicts the breaking scale for Gdark near ~ GeV. Alternatively, no fundamental scalar is needed, 
with the vector boson possibly generated by the condensate of a strongly coupled theory. 

At this juncture it is worth discussing an important point. As we have emphasized, to produce a Sommerfeld 
enhancement, we need an attractive interaction. Scalars (like gravitons) universally mediate attractive forces, but 
gauge fields can give attraction or repulsion, so do we generically get a Sommerfeld enhancement? For the case of 
the Majorana fermion (or real scalar) in particular, the dark matter does not carry any charge, and there does not 
appear to be any long-range force to speak of, so the question of the enhancement is more interesting. As we discuss 
in a little more detail in the appendix, the point is that the gauge symmetry is broken. The breaking dominates 
the properties of the asymptotic states. For instance, the dark matter must be part of a multiplet with at least 
two states, since a spin-1 particle cannot have a coupling to a single neutral state. The gauge symmetry breaking 
leads to a mass splitting between the states, which dominates the long-distance behavior of the theory, determining 
which state is the lightest, and therefore able to survive to the present time and serve as initial states for collisions 
leading to annihilation. However, if the mass splitting between the states is small enough compared to the kinetic 
energy of the collision, the gauge-partner DM states will necessarily be active in the collision, and eventually at 
distances smaller than the gauge boson masses, the gauge breaking is a negligible effect. Since the asymptotic states 
are in general roughly equal linear combinations of positive and negative charge gauge eigenstates, the remnant of 
asymptotic gauge breaking at short distances is that the incoming scattering states are linear combinations of gauge 
eigenstates, so the two-body wavefunction will be a linear combination of attractive and repulsive channels. While 
two-body wavefunction in the repulsive channel is indeed suppressed at the origin, the attractive part is enhanced. 
Therefore there is still a Sommerfeld enhancement, suppressed at most by an 0(1) factor reflecting the attractive 
component of the two-body state. (This is of course completely consistent with the Sommerfeld enhancements seen 




for ordinary WIMP annihilations, mediated by W/Z/7 exchange). 

Because of the presence of a new light state, the annihilation XX ~ * can > an d naturally will, be significant. In 
order not to spoil the success of nucleosynthesis, we cannot have very light new states in this sector, with a mass ^ 10 
MeV, in thermal equilibrium with the standard model; the simplest picture is therefore that all the light states in the 
dark sector have a mass ~ GeV. Without any special symmetries, there is no reason for any of these particles to be 
exactly stable, and the lightest ones can therefore only decay back to standard model states, indeed many SM states 
are also likely kinematically inaccessible, thus favoring ones that produce high energy positrons and electrons. This 
mechanism was first utilized in [l^ | to generate a large positron signal with smaller tt° and p signals. Consequently, an 
important question is the tendency of </> to decay to leptons. This is a simple matter of how couples to the standard 
model. (A more detailed discussion of this can be found in 30].) 

A scalar can couple with a dilaton-likc coupling ipF^ F^ Vl which will produce photons and hadrons (via gluons). 
Such a possibility will generally fail to produce a hard e + e~ spectrum. A more promising approach would be to mix 
with the standard model Higgs with a term K(p 2 h'h. Should acquire a vev (0) ~ m^, then we yield a small mixing 
with the standard model Higgs, and the will decay into the heaviest fermion pair available. For $ 200 MeV 
it will decay directly to e + e~, while for 200 MeV ^ m^ $ 250 MeV, will decay dominantly to muons. Above that 
hadronic states appear, and pion modes will dominate. Both e + e~ and fi + yT give good fits to the PAMELA data, 
while e+e~ gives a better fit to PAMELA+ATIC. 

A pseudoscalar, while not yielding a Sommcrfcld enhancement, could naturally be present in this new sector. Such 
a particle would typically couple to the heaviest particle available, or through the axion analog of the dilaton coupling 
above. Consequently, the decays of a pseudoscalar would be similar to those of the scalar. 

A vector boson will naturally mix with electromagnetism via the operator F' v F fiU . This possibility was considered 
some time ago in 40|. Such an operator will cause a vector </> M to couple directly to charge. Thus, for m^,^ Im^ it 
will decay to e + e~, while for 2m M ^$ ^ 2777,^ it will decay equally to e + e~ and /i + pT . Above 2771^, it will decay 
40%e + e~, 40% and 20%7r + 7r~. At these masses, no direct decays into 7r°'s will occur because they are neutral 
and the hadrons arc the appropriate degrees of freedom. At higher masses, where quarks and QCD are the appropriate 
degrees of freedom, the will decay to quarks, producing a wider range of hadronic states, including 7r°'s, and, at 
suitably high masses ^ 2 GeV, antiprotons as well 66] . In addition to XDM [3] , some other important exam ples 



41| , mirror dark matter [42 



of theories under which dark matter interacts with new forces include WIMPless models 
and secluded dark matter [4^ . 

Note that, while these interactions between the sectors can be small, they are all large enough to keep the dark 
and standard model sectors in thermal equilibrium down to temperatures far beneath the dark matter mass, and (as 
mentioned in the previous section), we can naturally get the correct thermal relic abundance with a weak-scale dark 
matter mass and perturbative annihilation cross sections. Kinetic equilibrium in these models is naturally maintained 



down to the temperature 7cmb ~ 44 1 . 
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IV. A NON-ABELIAN G da rk: INTEGRAL, DIRECT DETECTION, AND DAMA 

Up to this point we have focused on a situation where there is a single force-carrying boson 4>, whether vector or 
scalar. Already, this can have significant phenomenological consequences. In mixing with the standard model Higgs 
boson, there is a nuclear recoil cross section mediated by <j>. With technically natural parameters as described in [18| , 
the rate is unobservable, although in a two- Higgs doublet model the cross section is within reach of future experiments 



45|. 



In contrast, an 800 GeV WIMP which interacts via a particle that couples to charge is strongly constrained. Because 
the <p boson is light and couples to the electromagnetic vector current, there are strong limits. The cross section per 



nucleon for such a particle is 43 1 



167rZ 2 a SM aDarke 2 Ainc 

CT = 7T-7 I 11 ) 



\ fl GeV N 



/ V137- 1 / V 938 MeV / 



32/ \A) VlO- 3 / V137- 1 / V938McV/ V m 4> 



x 1.8 x 10" 37 cm 2 , 



where aDark is the coupling of the <f> to the dark matter, e describes the kinetic mixing, /x ne is the reduced mass 
of the DM- nucleon system and asM is the standard model electromagnetic coupling constant. With the parameters 



above, such a scattering cross section is excluded by the present CDMS 46J and XENON 47J bounds by 6 orders 



of magnitude. However, this limit can be evaded by splitting the two Majorana components of the Dirac fermion 



[22j or by splitting the scalar and pseudoscalar components of a complex scalar [48|, |49j . Since the vector coupling- 
is off-diagonal between these states, the nuclear recoil can only occur if there is sufficient kinetic energy to do so. If 
the splitting 8 > v 2 fj,/2 (where fi is the reduced mass of the WIMP-nucleus system) no scattering will occur. Such 
a splitting can easily arise for a U(l) symmetry by a U(l) breaking operator such as jjip c iph* h* which generates a 



small Majorana mass and splits the two components (see 23j for a discussion) 



Remarkably, for 5 ~ 100 keV one can reconcile the DAMA annual modulation signature with the null results of 
other experiments 2^, IE 24], in the "inelastic dark matter" scenario. We find that the ingredients for such a scenario 



occur here quite naturally. However, the splitting here must be 0(100 keV) and the origin of this scale is unknown, a 
point we shall address shortly. 

Exciting Dark Matter from a Non-Abelian Symmetry 

One of the strongest motivations for a ~ GeV mass (f> particle prior to the present ATIC and PAMELA data was 
in the context of eXciting dark matter [3| . In this scenario, dark matter excitations could occur in the center of the 
galaxy via inelastic scattering \X ~ * XX* ■ If <5 = m-x* — m x ~ 2m e , the decay \* ~^ X e+e ~ can generate the excess of 



511 keV x-rays seen from the galactic center by the INTEGRAL [50j, l5l| satellite. However, a large (nearly geometric) 
cross section is needed to produce the large numbers of positrons observed in the galactic center, necessitating a boson 
with mass of the order of the momentum transfer, i.e., ^ M x v ~ GeV, precisely the same scale as we require for 
the Sommerfeld enhancement. But where does the scale S ~ MeV come from? Remarkably, it arises radiatively at 



precisely the appropriate scale [67 1. 



We need the dark matter to have an excited state, and we will further assume the dark matter transforms under 
a non-Abelian gauge symmetry. Although an excited state can be present with simply a U(l), this only mediates 
the process XX ~ * X*X* ■ If this requires energy greater than Am e it is very hard to generate enough positrons to 
explain the INTEGRAL signal. If we assume the dark matter is a Majorana fermion, then it must transform as a 
real representation of the gauge symmetry. For a non-Abelian symmetry, the smallest such representation will be 
three-dimensional [such as a triplet of SU(2)]. This will allow a scattering X1X1 ~~ > X2X3- If m 3 is split from m-i ~ m\ 
by an amount S ~ MeV, we have arrived at the setup for the XDM explanation of the INTEGRAL signal. 

Because the gauge symmetry is Higgsed, we should expect a splitting between different states in the dark matter 
multiplet. This could arise already at tree- level, if the dark matter has direct couplings to the Higgs fields breaking 
the gauge symmetry; these could naturally be as large as the dark gauge breaking scale ~ GeV itself, which would 
be highly undesirable, since we need these splittings to be not much larger than the DM kinetic energies in order to 
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get a Sommerfeld enhancement to begin with. However, such direct couplings to the Higgs could be absent or very 
small (indeed most of the Yukawa couplings in the standard model are very small). We will assume that such a direct 
coupling is absent or negligible. The gauge breaking in the gauge boson masses then leads, at one loop, to splittings 
between different dark matter states, analogous to the familiar splitting between charged and neutral components of 
a Higgsino or Dirac neutrino. These splittings arise from infrared effects, and so are completely calculable, with sizes 
generically 0(amz) in the standard model or 0(am ( j > ) ~ MeV in the case at hand. Thus we find that the MeV 
splittings needed for XDM arise automatically once the mass scale of the (j> has been set to 0(GeV). 

We would like \i an d Xi to stay similar in mass, which can occur if the breaking pattern approximately preserves 
a custodial symmetry. However, if they are too degenerate, we are forced to take e < 10~ 5 in order to escape direct- 
detection constraints. On the other hand, we do not want too large of a splitting between X2 and Xi> as this would 
suppress the rate of positron production for INTEGRAL. Thus, we are compelled to consider 621 ~ 100 — 200 keV, 
which puts us precisely in the range relevant for the inelastic dark matter explanation of DAMA. (See Fig. 2J) 

All of these issues require detailed model-building, which we defer to future work; however, existence proofs are 
easy to construct. The Lagrangian is of the form 

£ = -CsM + ^Dark + Anix- (12) 

As a familiar example, imagine that Goark = SU(2) x U(l), with gauge bosons w^i and 6 M which we collectively refer 
to as a^i, and the dark matter multiplet x transforming as a triplet under SU(2) and neutral under the f (1); it could 
be either a scalar or fermion. We also assume that some set of Higgses completely break the symmetry. Working in 
unitary gauge, the tree-level dark sector Lagrangian is 

^Dark = ^Gaugc Kin. + ^ '''f, <'<''<*< •' + ' ' ' (13) 

where the m^- makes all the dark spin-1 fields massive, and • • • refers to other fields such as the physical Higgses 
that could be present. At one loop, this broken gauge symmetry will induce splittings between the 3 real DM states 
all of O (a Darken Dark) as just discussed above. The leading interaction between the two sectors is via kinetic mixing 
between the new U(l) and the photon (which is inherited from such a mixing with hypercharge): 

£mix = \^ V F^ (14) 

We put an e in front of this coupling because it is natural for this coupling to be small; it can be induced at one loop by 
integrating out some heavy states (of any mass between the GeV and Planck scales) charged under both this new U(l) 
and hypercharge. This can easily make e ~ 10 -4 — 10~ 3 . Even without a U(l), a similar mixing could be achieved 
with an "S-parameter" type operator Tr[($/M YG^F^, where $ is a dark Higgs field with quantum numbers such 
that ($/M) p transforms as an adjoint under G*Dark; it is reasonable to imagine that the scale M suppressing this 
operator is near the weak scale. 

Going back to PAMELA/ATIC, the non-Abelian self-couplings of the vector bosons can have an interesting effect 
on the annihilation process. For large enough ctDark, the gauge bosons radiate other soft and collinear gauge bosons 
leading to a "shower"; this can happen when aoark ^>g (M x /m^) ~ 1- While for quite perturbative values of aDark 
this is not an important effect, it could be interesting for larger values, and would lead to a greater multiplicity of 
softer e + e _ pairs in the final state. 

Before closing this section, it is important to point out that it is not merely a numerical accident that the excited 
states relevant for INTEGRAL and DAMA can actually be excited in the DM-DM collisions and DM collisions 
with direct-detection nuclei, but is rather a parametric consequence of maximizing the Sommerfeld enhancement for 
the annihilation cross section needed to explain ATIC/PAMELA. As we already emphasized in our discussion of 
Sommerfeld enhancement with vector states, vector boson couplings necessarily connect different dark matter mass 
eigenstates, and therefore there is no enhancement for the annihilation cross section needed to explain ATIC/PAMELA 
if the mass splittings are much larger than the kinetic energy available in the collision. But this parametrically implies 
that dark matter collisions should also have the kinetic energy needed to create the excited states, as necessary for 
the INTEGRAL signal. It is also interesting to note that the condition needed for the large geometric capture cross 
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X3 [i^x + i am 4> ~ lMeV)] 



X2 [m x + (a Try ~ lOOkeV)] 
Xi [m x ] 

FIG. 4: Spectrum of exciting dark matter. 

section, ^ Mv, also tells us that the Sommerfeld enhancement is as large as can be at these velocities, and has not 
yet been saturated by the finite range of the force carrier. Furthermore, since the mass of the heavy nuclei in direct 
detection experiments is comparable to the dark matter masses, the WIMP-nucleus kinetic energy is also naturally 
comparable to the excited state splittings. In this sense, even absent the direct experimental hints, signals like those 
of INTEGRAL and inelastic scattering for direct detection of dark matter are parametric predictions of our picture. 




V. SUBSTRUCTURE AND THE SOMMERFELD ENHANCEMENT 

As we have seen, the Sommerfeld enhancement leads to a cross section that scales at low energies as crv ~ 1/v. 
This results in a relatively higher contribution to the dark matter annihilation from low velocity particles. While 
the largest part of our halo is composed of dark matter particles with an approximately thermal distribution, there 
are subhalos with comparable or higher densities. Because these structures generally have lower velocity dispersions 
than the approximately thermal bulk of the halo, the Sommerfeld-enhanced cross sections can make these components 
especially important. 

Subhalos of the Milky Way halo are of particular interest, and N-body simulations predict that many should be 
present. There is still some debate as to what effect substructures can have on indirect detection prospects. However, 



some of these subhalos will have velocity dispersions of order 10 km/s 52|, and a simple examination of Fig. [5] shows 
that dramatic increases of up to two orders of magnitude in the Sommerfeld enhancement can occur for these lower 
velocities! 

Although in our model there are no direct ir° gammas, there are still significant hopes for detection of dwarf galaxies 
through dark matter annihilation. If the DM is also charged under SU(2) x U(l) there may be a subdominant 
component of annihilations into W + W~ for instance, much larger than the s-wave limited thermal cross section, 



which could yield significant photon signals from the hadronic shower 68j . The copious high energy positrons and 
electrons, produced even more abundantly than expected for a nonthermal WIMP generating the ATIC or PAMELA 
signal, can produce inverse Compton scattering signals off of the CMB. Even a loop suppressed annihilation into 77 
may be detectable with this enhancement. 

This presumes that the local annihilation does not have significant enhancement from a low- velocity component, 
however. If the local annihilation rate is also enhanced by substructure, then our expectations of the enhancement 
for dwarf galaxies would not be as large. As a concrete example, we can consider the possibility of a "dark disk." 
Recently, it has been argued that the old stars in the thick Milky Way disk should have an associated dark matter 
component, with dynamics which mirror those stars, and a density comparable to the local density [5^]. Because of 
the low velocity dispersion (of order 10 km/s), our estimates of what is a reasonable Sommerfeld boost may be off by 
an order of magnitude. 

As a consequence of this and other substructure, the Sommerfeld boost of Fig. [1] should be taken as a lower bound, 
with contributions from substructure likely increasing the cross section significantly further. We will not attempt to 
quantify these effects here beyond noting the possible increases of an order of magnitude or more suggested by Fig. 
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FIG. 5: Contours for the Sommerfeld enhancement factor S as a function of the mass ratio m,j > /m x and the coupling constant 
A, a = 10 km/s. 



[5j However, understanding them may be essential for understanding the size and spectrum if this enhancement of the 
cross section is responsible for the signals we observe. 



VI. OUTLOOK: IMPLICATIONS FOR PAMELA, PLANCK, FERMI/GLAST, AND THE LHC 

If the excesses in positrons and electrons seen by PAMELA and ATIC are arising from dark matter, there are 
important implications for a wide variety of experiments. It appears at this point that a simple modification of a 
standard candidate such as a neutralino in the minimal supersymmetric standard model (MSSM) is insufficient. The 
need for dominantly leptonic annihilation modes with large cross sections significantly changes our intuition for what 
we might see, where, and at what level. There are a few clear consequences looking forward. 

• The positron fraction seen by PAMELA should continue to rise up to the highest energies available to them 
(~ 270 GeV), and the electron + positron signal should deviate from a simple power law, as seen by ATIC. 

• If the PAMELA/ATIC signal came from a local source, we would not expect additional anomalous electronic 
activity elsewhere in the galaxy. Dark matter, on the other hand, should produce a significant signal in the 
center of the galaxy as well, yielding significant signals in the microwave range through synchrotron radiation 
and in gamma rays through inverse-Compton scattering. The former may already have been seen at WMAP 

h], HI and the latter at EGRET [H 1 141 ] . Additional data from Planck and Fermi/GLAST will make these 



signals robust [54 1 



• We have argued that the most natural way to generate such a large signal is through the presence of a new, light 
state, which decays dominantly to leptons. It is likely these states could naturally be produced at the LHC in 
some cascade, leading to highly boosted pairs of leptons as a generic signature of this scenario ■ 



• Although the cross section is not Sommerfeld enhanced during freeze-out, it can keep pace with the expansion 
rate over large periods of the cosmic history between kinetic decoupling and matter-radiation equality. This 
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can have significant implications for a variety of early-universe phenomena as well as the cosmic gamma-ray 
background Q. 

• The Sommerfeld enhancement is increasingly important at low velocity. Because substructure in the halo 
typically has velocity dispersions an order of magnitude smaller than the bulk of the halo, annihilations can be 
an order of magnitude higher, or more. With the already high local cross sections, this makes the prospects 
for detecting substructure even higher. With the mass range in question, continuum photons would be possibly 
visible at GLAST, while monochromatic photons [which could be generated in some models would be 
accessible to air Cerenkov telescopes, such as HESS (see 55 1 for a discussion, and [56| for HESS limits on DM 
annihilation in the Canis Major overdensity)]. 

We have argued that dark matter physics is far richer than usually thought, involving a multiplet of states and a 
new sector of dark forces. We have been led to propose this picture not by a flight of fancy but rather directly from 
experimental data. Even so, one can justifiably ask whether such extravagances are warranted. After all, experimental 
anomalies come and go, and it is entirely possible that the suite of hints that motivate our proposal are incorrect, or 
that they have more conventional explanations. However, we are very encouraged by the fact that the theory we have 
presented fits into a very reasonable picture of particle physics, is supported by overlapping pieces of experimental 
evidence, and that features of the theory motivated by one set of experimental anomalies automatically provide the 
ingredients to explain the others. Our focus in this paper has been on outlining this unified picture for dark matter; 
new experimental results coming soon should be able to tell us whether these ideas are even qualitatively on the right 
track. Needless to say it will then be important to find a specific and simple version of this theory, with a small 
number of parameters, to more quantitatively confront future data. 
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APPENDIX A: A QUICK REVIEW OF SOMMERFELD ENHANCEMENT 

The Sommerfeld enhancement is an elementary effect in nonrelativistic quantum mechanics; in this appendix we 
will review it in a simple way and discuss some of the parametrics for how the enhancement works for various kinds 
of interactions. 

Consider a nonrelativistic particle moving around some origin. There is an interaction Hamiltonian H ann = 
U ann S 3 (r) localized to the origin, which e.g. annihilates our particle or converts to another state in some way. 
Imagine that the particle is moving in the z direction so that its wavefunction is 

^°>(f) = e«* (Al) 

then the rate for this process is proportional to |V' < ' - ) (0)| 2 - But now suppose we also have a central potential V(r) 
attracting or repelling the particle to the origin. We could of course treat V perturbatively, but again at small 
velocities the potential may not be a small perturbation and can significantly distort the wave-function, which can be 
determined by solving the Schrodinger equation 

1 k 2 

V 2 ^fc + V(rU k = Vfe (A2) 

2M ' 2M y ' 
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with the boundary condition enforcing that the perturbation can only produce outgoing spherical waves as r — > oo: 

ikr 

ik-> 



ip -» e Ifcz + f{6) asr^oo (A3) 



r 



Now, since the annihilation is taking place locally near r = 0, the only effect of the perturbation V is to change the 
value of the modulus of the wave-function at the origin relative to its unperturbed value. Then, we can write 

CT = <T S k (A4) 

where the Sommerfeld enhancement factor S is simply 

Sk = = |^(0)|" (A5) 



where we are using the normalization of the wavefunction ipk as given by the asymptotic form of Eq. 

Finding a solution of the Schrodinger equation with these asymptotics is a completely elementary and standard 
part of scattering theory in nonrelativistic QM, which we quickly review for the sake of completeness. Any solution 
of the Schrodinger equation with rotational invariance around the z axis can be expanded as 

ij k =^AiPi{cos6)R H {r) (A6) 
i 

where Rki (r) are the continuum radial functions associated with angular momentum I satisfying 

1 1 d ( 2 d \ f 1(1 + 1) , A fc 2 /A s 

- ^7ZSl=l r 2 -R kl + ±±J. + V(r) Ru = r h (A7 ) 



2Mr 2 dr \ dr ' J \ r 2 v ' ) 2M 

The R k i(r) are real, and at infinity look like a spherical plane wave which we can choose to normalize as 

Rki(r) -» -sin(kr- -Ix + S^r)) (A8) 
r 2 

where 5i(r) <C kr as r — > oo. The phase shift Si(r) is determined by the requirement that Rki(r) is regular as r — > 0. 
Indeed, if the potential V(r) does not blow up faster than \/r near r — > 0, then we can ignore it relative to the kinetic 
terms, and we have that R k i(r) ~ r l as r — > 0; all but the / = terms vanish at the origin. We now have to choose 
the coefficients Ai in order to ensure that the asymptotics of Eq. [A3] are satisfied. Using the asymptotic expansion of 



gikz 



e kz ^ 1 y^, 2l + i\ Pl r cose \ \ e ikr _ (_!)J e -<fcr] (A9) 
2ikr £ — ' 



determines the expansion to be 

V^fc = lY. 1 ^ 21 + iy Sl Pi(cos9)R kl (r) (A10) 
l 

It is now very simple to determine ^(0), since as we just commented, R k i(r = 0) vanishes for every term other than 
1 = 0. Thus, we have 



Sk = 



Rk,i=o{0) 2 



(All) 



We can furthermore make the standard substitution Rk,i=o = Xk/ r > then the Schrodinger equation for \ turns into a 
one-dimensional problem 
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which we normalize at infinity with the condition 

Xfe(r) ->sin(fcr + <5) (A13) 
and since Rk,i=o goes to a constant as r — > 0, we have to have that x - > as r — > 0, or 

Xfc(r)-»r^(0) asr^O (A14) 

Effectively, we can imagine launching % from % = at r = with different velocities ^p-(O), and these will evolve to 
some waveform as r — > oo, but the correct x'(0) is determined by the requirement that the waveform at infinity have 
unit amplitude. 



Summarizing, then, the Sommerfeld enhancement is 



(A15) 



where Xfc satisfies the ID Schrodinger equation [XT2] with boundary conditions Eqs lA13llA141 As a sanity check, let us 
see how this works with vanishing potential. The solution that vanishes as r — > is Xfc( r ) — ^lsin(fcr), and matching 
the asymptotics forces A = 1. Then x' fc (0) = k and Sk = 1. 

In the recent literature on the subject, a different expression for the Sommerfeld enhancement is used, arising from 
the use of the optical theorem. We are instructed to solve the same ID Schrodinger equation (Eq. IA12[) . this time 
with no special boundary conditions at x = 0, but with boundary conditions so that x oc e +tkr as r — ► oo. Then, the 
Sommerfeld enhancement is said to be 

Sk = "boWF (A16) 

It is very easy to show that that these two forms for Sk agree exactly. To see this, let us begin by denoting xi( r ) to 
be the solution to the Schrodinger equation (Eq. IAl2j) with the boundary condition xi( r ) — > sin(fcr + 5) as r — > oo 
(Eq. IA13[) . As shown above, Xi(0) = 0. Let X2M be the linearly independent solution with the boundary condition 
X2( r ) — * cos(fcr + S) as r — > 00, and define A = X2(0). Now Eq. IA12l has a conserved Wronskian 

W = X i(r)X2(r)-X2(r)x'i(r). (A17) 

It is easy to verify directly from the differential equation that W'(r) = 0; this is true (Abel's theorem) because there 
are no x' terms in the differential equation. But comparing the values of the conserved Wronskian at zero and 00, 

W(oo) = -fc(sin 2 (fcr + S) + cos 2 (fcr + S)) = -k = W(0) = -^x'i(°)- ( A18 ) 

So then |xi(0)| = and our new expression for the Sommerfeld enhancement, Eq. IA151 is just Sk — 1/\A\ 2 . Now, 

our second form Sk = |x(°°)| 2 > where x( r ) satisfies the boundary conditions x'( r ) ~^ ^x{ r ) as r 00, and x(0) = E 
By the asymptotic behavior at large r, we can identify x( r ) as the linear combination x( r ) — C{X2( r ) +*Xi( r ))7 where 
C is some complex constant. But then x(0) = CA — 1, and Sk — \C\ 2 , so as previously we obtain Sk = 1/\A\ 2 . Thus 
the two formulae for the Sommerfeld enhancement are equivalent. 



1. Attractive Coulomb Potential 

Let us see how this works in some simple examples. We are solving the Schrodinger equation for a particle of mass 
M and asymptotic velocity v, with potential 

V(r) = -|- (A19) 
which we solve with the boundary conditions of Eqs. IA13IA141 
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We simplify the analysis by working with the natural dimensionless variables, with the unit of length normalized 
to the Bohr radius, i.e. we take 



Then the Schrodingcr equation becomes 



where we have defined the parameter 



r = cT x M- x x (A20) 



- x" ~\x = £ 2 X (A21) 



e v = - (A22) 
a 



Of course we can solve this exactly in terms of hypergeometric functions to find the result obtained by Sommerfcld 



Sk = 



1 - e 



(A23) 



Note that as e v — > oo, Sk — » 1 as expected; there is no enhancement at large velocity. For the attractive Yukawa at 
small velocities we have the enhancement 

S k - — (A24) 

v 

while for the repulsive case, there is instead the expected exponential suppression from the need to tunnel through 
the Coulomb barrier 

S k ~ e'^ (A25) 

To get some simple insight into what is going on, let us re-derive these results approximately. For x much smaller 
than 1/e 2 , we can ignore the kinetic term. In the WKB approximation, the waves are of the form x 1 / 4 e lv/ ^, so the 
amplitudes grow like x 1 ^ 4 . In order to match to a unit norm wave near x ~ we have to scale the wavefunction 

1/2 

at small a; by a factor ~ e v ' , so that near the origin 

X ~ xe x J 2 ~ aMrel' 2 (A26) 
from which we can read off the derivative at the origin &(0) = el /2 aM, and with k = Mv we can determine 

^H^! 2 = - (A27) 
Mv v 

which is correct parametrically. We can also arrive at this result from the second form for Sk- The computation is 
even more direct here. We have waveforms growing as x 1 / 4 towards x ~ 4r- In the region near x <~ l/e v , we must 
transition to a purely outgoing wave. This is a transmission/reflection problem, and ingoing and outgoing waves from 
the left will have comparable amplitude. When we continue these back to the origin, we will then have an amplitude 
reduced by a factor ~ ej 2 . Then, 

S k ^(4j,) 2 -^ (A28) 



2. Attractive Well and Resonance Scattering 

Let us do another example, where 



V = -V o 0(L-r),V Q =^ (A29) 
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and so 



The solution inside is x(r) = Asinki n r, where kf n = n 2 + k 2 , while outside we write it as sin(fc(r — L) + S. Then, 
matching across the boundary at r = L gives 

A sin ki n L = sin S, k in A cos k in L — k cos 5 (A30) 

Squaring these equations and adding them we can determine 

A 2 = — \r (A31) 

sin ki n L H — p 2 - cos ki n L 

A 2 k 2 1 
k A- sin ki n L + cos z ki n L 

Now for k 2 /2M <C Vo, we have k in — k + /c 2 /(2k) + - • • . Clearly, if coskL is not close to zero, there is no enhancement. 
However, if cos nL = 0, then we have a large enhancement 

Sk - ^ (A33) 

This has a very simple physical interpretation in terms of resonance with a zero-energy bound state. Our well has a 
number of bound states, and typically the binding energies are of order Vo- We see that if coskL is not close to zero, 
we have to have A <~ k/ki n be small. However, if accidentally cos kL = 0, then there is a zero-energy bound state: 
the wavefunction can match on to ip = 1 for r > L smoothly, giving a zero-energy bound state. The enhancement is 
of the form 

S - g (A34) 

E - Abound k 2 

Note this formally diverges as v — > 0, but is actually cut off by the finite width of the state as familiar from Breit- 
Wigner. 



Now let us examine 



Working again in Bohr units, we have 



Attractive Yukawa Potential 



a 



V(r) = e- m * r 



(A35) 



V(x) 



(A36) 



where 



m 
~ aM' 



(A37) 



Now, if -c e 2 , the Yukawa term is always irrelevant and we revert to our previous Coulomb analysis. 

However, if » e 2 , our analysis changes; we will use the second expression for the Sommerfeld enhancement for 
simplicity. The potential turns off exponentially around x <~ l/e^. Now, the effective momentum is 



-e-^ x + el 



(A38) 



and the quantity 



ft cff 
k 2 



(A39) 
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determines the length scale the potential is varying over relative to the wavelength; so long as it is small, the WKB 
approximation is good, and we have a waveform growing as k e ^ 2 e 1 ^ dx feef f( x ). Note that for 1 -C x <C l/e<£, the 
WKB approximation is manifestly good. Let us now take the arbitrarily low velocity limit, where e v — > 0. Then in 
the neighborhood of x ~ l/e^ we have ~ e^e -6 *^, and 

fete***' - 60 



ft cff 



(A40) 

1/2 

' . The potential then 

varies more sharply than the wavelength, and we have a reflection/transmission problem, with an 0(1) fraction of the 
amplitude escaping to infinity. The enhancement is then 



so the WKB approximation breaks down when k c g ~ e^, where the WKB amplitude is ~ e 



aM 



(A41) 



We did this analysis for e v — > 0, but clearly it will hold for larger e„, till e v ~ e^, at which point it matches smoothly 
to the j- enhancement we get for the Coulomb problem. The crossover with e v ~ is equivalent to Mv ~ m$, when 
the deBroglie wavelength of the particle is comparable to the range of the interaction. This is intuitive-as the particle 
velocity drops and the deBroglie wavelength becomes larger than the range of the attractive force, the enhancement 
saturates. Of course if is close to the values that make the Yukawa potential have zero-energy bound states, then 
the enhancement is much larger; we can get an additional enhancement ~ e^/e^ up to the point where it gets cut off 
by finite width effects. 

In this simple theory it is of course also straightforward to solve for the Sommerfeld enhancement numerically. We 
show the enhancement as a function of and e v in Figs. [6] and [3 
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FIG. 6: Contour plot of S as a function of and e v . The lower right triangle corresponds to the zero-mass limit, whereas the 
upper left triangle is the resonance region. 



4. Two-particle annihilation 



Let us finally consider our real case of interest, involving two-particle annihilation. To keep things simple, let us 
imagine that the two particles are not identical, for instance they could be Majorana fermions with opposite spins; we 
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FIG. 7: As in Fig. [6] showing the resonance region in more detail. 



can restrict to this case because none of the interactions we consider depend on spin, and the annihilation channels 
we imagine can all proceed from this spin configuration. Let us imagine that there are a number of states \A of 
(nearly) equal mass, with the label A running from A = 1, - • • ,N. The two-particle Hilbert space is labeled by 
states \x, A;y, B) and a general wavefunction is {ip\x,A;y,B) = ipAB(x,y). There is some short-range annihilation 
Hamiltonian U6 3 (xi — x%) into decay products \ final); where U is an operator 



(final\U\AB) = M 



an a 
AB 



(A42) 



Now, suppose there is also a long-range interaction between the two particles with a potential. The general Schrodinger 
equation is of the form 



1 

2M 



(Va + Vy) ipAB{x,y) + V A bcd(x - y)ipcD(x,y) = Ei(j AB (x,y) 



As usual we factor out the center-of-mass motion by writing i^ABix^y) = e lP ' i, - xJr ^ ())ab(x — y), and we have 

1 



A I 



V r 4>AB(r) + VABCD(r)(t)CD(r) = EcM<j>AB{r) 



(A43) 



(A44) 



where Ecu is the energy in the center-of-mass frame. 

We have an initial state with some definite A* , , and we take as an unperturbed solution 

< j,^B')(p) AB = 5 A. s B' e ik, 

The annihilation cross section without the interaction V is proportional to 
However with the new interaction, the annihilation cross section is 



(A45) 



(A46) 



<JannK\^ B '\0)M CD \ 2 



(A47) 
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So we can write 

<W = x Sf B * (A48) 

where 

|*&f(oo)M£Sf [ yj 

Just as above, in reducing the problem to an S'-wave computation, we can replace 4>c® (0) with (x'( r — ®))cd m 
the obvious way. 

Of course one contribution to Vabcd comes from the small mass splittings between the states. If we write the 
(almost equal) common mass term for the DM states as (M + AMa)xaXa, the AM's show up in the potential as 

V s a P bcd = (AMa + AM B )S AC S BD (A50) 

For the Sommerfeld enhancement, we also need some long-range attractive interaction. As we have discussed, vectors 
are possibly the most promising candidate. The leading coupling to spin-1 particles a^i is 

9XA^XBT AB a^ (A51) 

and ignoring the mass of the gauge boson this gives us a 1/r contribution to the effective potential 

VTbTd (r) = -a-T AC T BD (A52) 

while taking the vector masses into account gives both Yukawa exponential factors and a more complicated tensor 
structure. In total, 

Vabcd = V*« D + vr£ e D (A53) 
Now, it is obvious that in our basis, the T AB should be antisymmetric 

Tab = -T BA (A54) 

since the gauge symmetry must be a subgroup of the SO(N) global symmetry preserved by the large common mass 
term M\aXa, and the SO(N) generators are antisymmetric. Thus, the coupling to vectors in this basis is necessarily 
off-diagonal. Let us look at a simple example, where N = 2 and we have a single Abelian gauge field. The 4>ab span a 
4-dimensional Hilbert space, with |11), |12), |21), 1 22} as a basis. Since the gauge boson exchange necessarily changes 
1 <-> 2, Vabcd is block diagonal, operating in two separate Hilbert spaces, spanned by (|11), 1 22} ) and ( 1 1 2} , |21)). 
Since we are ultimately interested in scattering with 11 initial states, let us look at the first one, where we have 

in the basis |11) = ^ ^ , |22) = ^ q ) ' ^ A55 ^ 

Clearly, if AM is enormous, we will not have any interesting Sommerfeld enhancement in (11) scattering, since in this 
case there is no long-range interaction between 11 at all, so let us assume that AM is smaller than the kinetic energy 
of the collision. Now it is clear that as r — > oo, the mass splitting dominates the potential, and obviously particle 1 is 
the lightest state. However, at smaller distances, the gauge exchange term dominates. This is not diagonal in the same 
basis, and has "attractive" and "repulsive" eigenstates with energies ±a/r. Note that the asymptotic |11) state is an 
equal linear combination of the attractive and repulsive channels. While the repulsive channels suffer a Sommerfeld 
suppression, the attractive channel is Sommerfeld enhanced. Note that as long as AM is paramctrically smaller than 
the kinetic energy, its only role in this discussion was to split the two asymptotic states, and therefore determine what 
the natural initial states are. Note also that if AM is large but not infinite, the mixing with 2 generates an attractive 
potential between 11 of the form V e ff(r) = —a 2 /(AMr). It would be interesting to understand these limits of the 
multistate Sommerfeld effect in parametric detail; we defer this to future work. 
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It is very easy to see that our conclusion about the presence of a Sommcrfeld effect is general for any gauge 
interaction. Wc think of V(ab)(cd) as a ma tri x hi the Hilbert space spanned by (AB). Note that since the T % 
are antisymmetric, they are also traceless, and as a consequence, the matrix V^b)(cd) * s a ^ so Graceless and so has 
both positive and negative eigenvalues, reflecting the obvious fact that gauge exchange gives us both attractive and 
repulsive potentials. Let us go to a basis in (AB) space where V 9auge is diagonal, and denote eigenvectors with negative 
(attractive) eigenvalues as f attractive an( j re p U l s i ve ones as j^v ulsme _ ^ e can think f t ne initial wavefunction in AB 
space <5^ B B as a state in (AB) space and expand it in terms of these eigenvectors as 

c(A*B*) _ ^A" B* f attractive , /~ iA * B * f repulsive / \ ro\ 

°AB — u ot( JAB 1" Wep JAB {AOOj 

wc can determine the coefheients by dotting left-hand side and right-hand side into the eigenvectors, so that 

c(A*B") f attractive f attractive , rrepulsive r repulsive /ir« 

°AB — JA»B* JAB ' J A* B* JAB \ A0 ' J 

Then, since the repulsive components are exponentially suppressed at the origin while the attractive components are 
enhanced, we get a Sommerfeld enhancement as long as /a* b* ^ 0. In particular, for scattering the same species, this 
is true so long as ^ 0. Said more colloquially, these Majorana states are linear combinations of "positive" and 

"negative" charged states; so long as they have any component which would mutually attract, there is a Sommerfeld 
enhancement from that component alone. 
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Although, it should be noted that such an explanation seems to fail because the spectrum would not be sufficiently hard; 
see Fig. 6 of [ll[. 

The importance of this effect in the context of multi-TeV scale dark matter interacting via W and Z bosons was first 
discussed by (j^ . and more recently emphasized with regards to PAMELA in the context of "minimal dark matter," with 
similar masses by [§3, [5^] . 

Another possibility for sufficiently light bosons would be radiative capture [Hit ]. 

Recently, [do| explored the possibility that PAMELA could distinguish the spin of the dark matter, which could be relevant 
here in annihilations to vector bosons. 
[67] Note that in the early universe, all these states will be populated at freeze-out, and will participate in maintaining thermal 
equilibrium. However, unlike coannihilation in the MSSM, which typically occurs between a Bino and a slepton, there is 
not expected to be a large effect, on the connection between the present annihilation rate and that at freeze-out. This 
is because the self-annihilation and coannihilation cross sections are similar or even equal here, while in the MSSM the 
coannihilation cross section is much larger. 
[68] Although larger annihilation rates into W + W~ or other standard model states are possible, the cross section must be 
lower than the limits of 0, [6l| arising from antiprotons. Moreover, strong limits from final state radiation [fj^] can also 
arise, although these depend more sensitively on the nature of the halo profile. 



